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We study nonperturbative multiphoton electron-positron pair creation in ultrastrong electromagnetic fields formed 
by two counterpropagating pulses with elliptic polarization. Our numerical approach allows us to take into account 
the temporal as well as the spatial variation of the standing electromagnetic field. The spin and momentum 
resolved pair creation probabilities feature characteristic Rabi oscillations and resonance spectra. Therefore, each 
laser frequency features a specific momentum distribution of the created particles. We find that, depending on 
the relative polarization of both pulses, the created electrons may be spin polarized along the direction of field 
propagation. 
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1. Introduction 

Quantum electrodynamics predicts the possible breakdown of 
the vacuum in the presence of ultrastrong electromagnetic fields 
into pairs of electrons and positrons. In the seminal articles 
by Sauter and others [1-3], pair creation was investigated the¬ 
oretically for constant electromagnetic fields. Since then this 
process has been studied in numerous works; see Refs. [4, 5] 
for recent reviews. The Schwinger critical field strength of 
Es = 1.3 X 10^^ V/m, where spontaneous pair creation is ex¬ 
pected to set in, cannot be reached even by the strongest laser 
facilities available today. However, pair creation may be as¬ 
sisted by additional fields or particles or by electromagnetic 
fields that oscillate in time and space. Novel light sources 
envisage to provide field intensities in excess of 10^^ W/cm^ 
and field frequencies in the x-ray domain [6-10]. The ELI- 
Ultra High Field Facility aims to reach intensities exceeding 
even 10^^ W/cm^ corresponding to a field strength of about 
10^^ V/m, which is only a few orders of magnitude below the 
critical field strength Es . It has been argued that the intensity 
of such near-future ultrastrong light sources may be sufficient 
to observe pair creation [11-15]. 

Pair creation in laser fields is often studied by considering a 
time-varying homogeneous electric field [12, 13, 16-22]. Ne¬ 
glecting the magnetic field component and the spatial variation 
of the electric field eases the theoretical analysis and is moti¬ 
vated by the fact that a standing electromagnetic wave, which 
is formed by superimposing two counterpropagating linearly 
polarized waves of equal wavelength, intensity, and polariza¬ 
tion, has a vanishing magnetic field component at positions, 
where the electric field reaches its maximum. A more realistic 
description of electron-positron pair creation, however, requires 
the consideration of inhomogeneous electric fields [23-27] or 
electromagnetic fields that depend on time and space [28-30]. 

Introducing the adiabaticity parameter^ = eEKmocoj) [21], 
where mo denotes the electron mass, c the speed of light, e the 
elementary charge, E the peak electric field strength, and oj the 


electromagnetic field’s angular frequency, one can distinguish 
two limiting cases of pair creation in oscillating electromagnetic 
fields [16] with frequencies oj < moc^lh. These two regimes 
are commonly referred to as nonperturbative Schwinger pair 
creation » 1) and perturbative multiphoton pair creation 
1). Note that these regimes are sometimes distinguished 
via the Keldysh parameter of vacuum pair creation y, which 
is the inverse of In this contribution we will focus on the 
intermediate regime of nonperturbative multiphoton pair cre¬ 
ation ^ 1), which is most challenging from the theoretical 
point of view and, therefore, less well understood [21, 28, 31]. 
In particular, we consider pair creation in a standing electro¬ 
magnetic wave formed by two counterpropagating waves. This 
means we include not only the electric field component but 
also the field’s magnetic component, which can alter the pair 
creation rate [32, 33]. The temporal as well as spatial oscilla¬ 
tions of these fields are accounted for. In our work the waves 
may have arbitrary elliptical polarization in contrast to prior 
publications where mainly linear polarization was considered. 
Electron dynamics in intense light with elliptical polarization 
[34] and pair creation via the nonlinear Bethe-Heitler process 
[35, 36] can show nontrivial spin effects [34, 35]. Therefore, 
pair creation in elliptically polarized light beams is expected to 
feature unexplored spin effects. 

The manuscript is organized as follows. In Sec. 2 we de¬ 
scribe the laser field configuration and pay special attention to 
the polarization of the light. We lay down the theoretical foun¬ 
dation for the numerical work in Sec. 3. The numerical methods 
are explained in Sec. 4. Our results are presented in Sec. 5. 
Typical Rabi oscillations of the pair creation probability and 
properties of the pair creation resonances are discussed. These 
resonance spectra are investigated for different polarizations 
and compared to a model that employs the dipole approxima¬ 
tion, i. e., neglecting the magnetic field component and the 
spatial variation of the electric field. The dependence of pair 
creation on the final momenta of the created particles is dis¬ 
cussed at the end of Sec. 5. Finally, we conclude in Sec. 6. For 
the reminder of this article, we will employ units with fi = I 
and c = 1. 
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2. Elliptically polarized light beams 

In this work we will focus on a setup with two counterpropagat- 
ing elliptically polarized laser beams with equal wavelength A. 
Note that two counterpropagating laser beams with different 
wavelengths can always be transformed by a suitable Lorentz 
boost along the propagation direction into the considered setup. 
Furthermore, two crossed laser beams with the same wave¬ 
length and arbitrary but nonzero collision angle can always be 
transformed into the considered setup by a suitable Lorentz 
boost perpendicularly to the propagation direction. Similarly, 
two crossed laser beams with different wavelengths and arbi¬ 
trary but nonzero collision angle can always be transformed by 
a suitable proper Lorentz transformation (which may include 
boosts as well as rotations) into the considered setup. Thus, it 
is quite generic. 

Let us denote the unit vectors along the coordinate axes of 
the Cartesian coordinate system by ejc, e^, and e^. We assume 
that the two light beams propagate along the z direction, and 
thus their wave vectors equal k+ = ±ke^ with k = InjA. The 
indices distinguish beams traveling to the right (index “-r”) and 
to the left (index respectively. To describe the polarization 

of an electromagnetic wave, it is convenient to use the Jones 
vector formalism [37]. The so-called Jones vectors |/) and \r) 
for circular polarized beams with left-handed and right-handed 
orientation are defined as 


\l)=A-{e^ + iey), (la) 

V2 

\r}=A-(e^-iey). (lb) 

V2 

Any polarization state of the electric fields of the two laser 
beams can be formed by a linear combination of the Jones 
vectors |/) and |r). At position r = (x,y,z) and time L the 
electric field of a single beam is given by 

E+(r, t) = Re (e (cos a+ \l) + sin a+ |r)) Q<±kz-cot)"^ ^ ^2) 

where co = k denotes the angular frequency and E the peak 
electric field strength. The expression (2) may be transformed 
into 


E±{r, t) = E 


Aos if+jl 
sin ^+/2 


- sin ^+/2' 
cos ^+/2 


X 


0 


/ 


/cos(+kz - cot + (fi±l2) cos(a+ - 7r/4)' 
y sin(+kz - cot (fi±l2) sin(Q(+ - 7r/4) 


( 3 ) 


The second factor in (3) parametrizes an ellipse, and the pa¬ 
rameters a+ and a- determine its ellipticity. The parameters 
(p+ and (f- specify the orientation of this ellipse via the first 
factor in (3), which represents a rotation matrix; see also Fig. 1. 
The light beams have linear polarization for a± = njA and 
a+ = 37:14, while they have circular polarization for q(+ = 0 
and also for a+ = 7r/2. The magnetic field component follows 
via B+ = k+x EJk, and the mean intensity equals soE^/2 in¬ 
dependent of a± and . The corresponding Coulomb-gauge 


a = 0 

»=o(3 




^ = 0 (f = t: 

^ 0 


FIG. 1 : (color online) Schematic illustration of different polarization 
states as a function of the parameters a and (p, which determine the 
ellipticity (degree of polarization) and the orientation of the polariza¬ 
tion ellipsis in the x-y plane. Upper row for constant orientation cp but 
varying a; lower row for constant ellipticity a but varying ip. 


a+ - a- rr+ = ^/2 - a- 


FIG. 2: (color online) Schematic illustration of the possible relative 
spin orientations of two counterpropagating electromagnetic fields. 
Left: When the fields’ ellipticity parameters equal (u_ = a+), then 
both fields have the same sense of rotation. Consequently, also the 
superposition of both fields rotates around the z axis and this axis of 
rotation leads to a preferred direction along the positive or negative z 
axis. Right: For a- = 7:12-a + the sense of rotation is exaetly opposite. 
Thus, the resulting field does not rotate, and no direction of the z axis 
is preferred. 


vector potential is 


^±iT, 0 = Re (cos a+ \l) + sin a+ e^*^- |r)) j (4) 

or equivalently 


A±{r,t) 


E 

CO 


'cos ip+12 
sin ipJ2 


- sin if ±12^ 
cos ip+12 


X 
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i sin(+kz - cot + ip+12) cos(a+ - 7r/4) 
cos(+kz - cot + (p+12) sin(o(+ - 7r/4) 


( 5 ) 


The total vector potential of two counterpropagating electro¬ 
magnetic waves equals A = A+ -r A_ and similarly for the 
electric and magnetic field components. In the reminder of 
this article, we will assume that both counterpropagating light 
waves have the same orientation, i. q., ip+ = cp-, and that they 
have the same ellipticity, Lq., a- = a+ or a- = 7r/2 - a+ de¬ 
pending on the sense of rotation being identical or opposite to 
each other; see Fig. 2. Note that for any such field configuration 
we can always find a coordinate system by a rotation around the 
z axis, where ip+ = ip- = 0. Thus, we will set in the following 
ip+ = ip_ = 0 without loss of generality. 

Elliptically polarized light beams carry photonic spin angu¬ 
lar momentum [38], which is sqE^ cos(2a+)l(2co) for each of 
the two considered counterpropagating light beams. Thus, the 
ellipticity paramaters a+ and a- determine the photonic spin 
density. The total spin density of the considered setup equals 
2£oE^ cos(2a+)lco (with orientation along the z direction) for 
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the case that both individual beams have the opposite helicity, 
i. e., a+ = a-; see also the left part of Fig. 2. The total photonic 
spin density vanishes if both light beams have same helicity, 
i. e., a+ = 71 12 - a-; see also the right part of Fig. 2. In what 
follows, we will name the field setup “corotating” if a+ = cr_ 
and “antirotating” if a+ = ti/I - a-. 


3. Theoretical foundations 

Quantum electrodynamics has to be applied to describe pair 
creation in electromagnetic fields [39]. The Hamiltonian de¬ 
scribing the dynamics of a quantum field state in an external 
electromagnetic field is given by 

H(t) = ^ 4f(rfm)Hr)dh, (6) 

where \J/{r) denotes the time-independent spinor field operator 
and 

7T(0 = or • ( - iV - qA{r, t)) + mofi (7) 

is the Dirac Hamiltonian for a particle with charge q and rest 
mass mo. Furthermore, a = and (3 stand for the 

Dirac matrices [40, 41]. Let us denote the times, when the 
interaction with the external electromagnetic field sets in and 
ends, by ^in and ^out, respectively. Thus, E(t) = 0 for r < ^in, 
and for t > ^out- The Dirac Hamiltonian may be diagonalized 
at times and ^out by introducing the two complete sets of 
orthonormal functions +^„(r) and -(fnix) [42], 

= ±Sn , (8a) 

'H(rout)“^n(^) = ~<Pnir) , (8b) 

where +£„ > 0, _£„ < 0 and > 0, “£„ < 0 denote the corre¬ 
sponding positive and negative eigenenergies, with n labeling 
the quantum state. Then the spinor field operator \j/{r) may be 
decomposed in one of the two bases 

Hr) = L +^n{r)an{tm) + -iPn{r)bl{tia ), (9a) 

n 

0^(^) ~ ^ ^n(^)^«(^out) "I" (9b) 

n 

where , d and denote the creation and the annihilation 
operators for the electron and the positron, respectively, at times 
tin and ^out- Quantum field states before and after interaction 
with definite numbers of particles can always be written as 
products of the vacuum states |0, ^in) and |0, ^out) and a number 
of particle and antiparticle creation operators, 

Ifin) — d^^^{tiTn)d\^{tin) . . . • • • |0? fin) 9 (10a) 

kout) ~ (^out)^l2^^out) • • • (^out)^m2^^out) • • • IO9 ^out) • 

(10b) 

The amplitude for the transition between these states is given 
by 

4 ^in—>out — (Lutl^(Lut9 fin)lfin) ? (H) 


where U fin) denotes the time evolution operator that corre¬ 
sponds to the Hamiltonian H{t) and maps the state I fin) at time 
fin to some later time ^out- 

After interaction the average number of electrons in quantum 
state n is given by 


K 

— {U (^out? fin)finl^^(^out)^n(^out)|f^(^out? fin)fin) 

(12) 

or equivalently 



— (finl<5^(Lut)<5n(Lut)lfin) ? 

(13) 

where 

dn(jout) ~ G (Lut? fin)^^«(^out)U(^out? fin) 

(14) 


is the representation of the electron creation operator in the 
Heisenberg picture with l^in) as the reference state. In the 
Heisenberg picture, the spinor field operator \J/{r) is replaced 
by the time-dependent operator ij/{r, t) of which the evolution 
is determined by the Dirac equation 

= (15) 

At time fin the operator t) equals (9a), and at later times 
Hr,l) = U{l,tiJHr)U{l,ti„) (16) 

or equivalently 

Hr, ^) = J' tonur', tin)Hr) dV', (17) 

where G(rJ\r\t') denotes the complete Green function of 
(15). Using 

4(?out) = ^ ^iPn(r)^(r,t)A^r, (18) 

(9a), and (17), we get 

dni3o\x\) — ^ f^( \+)n-,m^mihn) + f^( l-)/i;m^m(fin) (19) 

m 

with G{-\+)n-m defined as 

= *^„(r)^G(/-,rout;/-',6n)±^m(r')dV'dV. (20) 

Choosing the vacuum state |0, fin) as the initial state Ifin)? one 
can show with (19) that [43] 

Nn = Y, GC\-):,^GC\-)n-.m = ^ \GC\-)n-.mf , ( 21 ) 

m m 

and therefore the expectation value of the total number of cre¬ 
ated electrons is 

N=Y,\GC\-)n;J. ( 22 ) 

m,n 

Note that the quantity \GC\-)n-,m\^ denotes the probability that 
a negative-energy state with the quantum number m turns into 
a positive-energy state with the quantum number n. A similar 
calculation yields the expectation value of the total number 
of created positrons 2m,n \G(~\+)n-m\^- In order to keep the 
manuscript compact, we will present results for the created 
electrons only. Due to the synnnetry of the pair creation process, 
corresponding results for the positron can be deduced from 
calculations for the electron. 
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4. Numerical solution of the Dirac 
equation in monochromatic fields 

In Sec. 3 we outlined that the number of created pairs can 
be determined via the quantities G(-|+)„;^, which can also be 
written as 

Gt\±)n-m = J'''<Pn(r)K‘fim(r,tout)d^f, (23) 

where t) is a solution of the time-dependent Dirac equa¬ 
tion (15) with the initial condition ±(fm(r, ^m) = Thus, 

we can analyze the quantum electrodynamical problem of 
pair creation by solving a large number of independent time- 
dependent single-particle problems for the Dirac equation (15). 
For the setup of Sec. 2, it is appropriate to solve these in momen¬ 
tum space. The (real-valued) total vector potential A = A+ -r A_ 
can always be written in the form 

A(r, t) = (24) 

with some complex-valued vector a{t) and its conjugate a{ty. 
Because the external laser field was assumed to be monochro¬ 
matic, the operator a • A (and consequently also the Dirac 
Hamiltonian) couples momentum eigenstates with momentum 
p only to momentum eigenstates with momentum p±kez, which 
simplifies the numerical treatment considerably. Because of 
the monochromatic field (24), the Hilbert space separates into 
disjunct subspaces, and each of them is spanned by plane waves 
with momentum nke^ + po with integer n, where the base mo¬ 
mentum po selects a specific subspace. To take advantage of 
this symmetry, we expand the single-particle wave functions 
into plane waves as 


has been introduced to model a smooth turn on and turn off of 
the external laser field at ^ = tin = -AT and t = tout = T AT. 

As the vector potential vanishes at tin and ^out. the states 
±(Pn(r) and -^Pnix) are free-particle states with definite momen¬ 
tum, energy, and spin orientation (in the z direction), which 
read explicitly 

+V>n,vi{r) = ^(fnMitr) = , (30a) 

->PnMi{r) = >„.t/i(r) = , (30b) 

where we have replaced the generic quantum number n by the 
quantum numbers (n, t) or (n, J,), which specify the particle’s 
momentum {nke^ + po) and its spin state (up or down). We also 
introduced 


M, 


'p.vi - 


Vp,Vl = 


mo+6(p) 


I 2S(p) 
mo + S(p) 


26(p) 


Xvi 

(rp 

m + £>(p) 

o-p 

mo + Sip) 

Xvi 


'Xmi 


[Xvi 


(31a) 


(31b) 


with;i^l = (1,0)^ and;p| = (0,1)^ and 6(p) = P^‘ 

Fourier space the states +(pn,yi{r) and ^ipn,yi{r) correspond to 


di = 


I ^nke^+P(j,'\ll for I — 11 
(0,0,0,0)"^ else. 


while and tpnyuir) correspond to 

_ \^nke-+p(i,'[li for / — Tl 


di 


I (0,0,0,0)"^ else. 


(32a) 


(32b) 


«F(r, 0 = y dlit) J (25) 

^ y ZTT 

n,y 

with 7 € {1,2, 3,4} and the vectors 

-a'=( 1,0,0,0)T, (26a) 

= (0,1,0,0)"^, (26b) 

= (0,0,1,0)"^, (26c) 

lA"^ = (0,0,0, l)"^. (26d) 

The expansion (25) into a discrete set of plane waves would 
also be justified for multicolor setups provided that all wave 
vectors are an integer multiple of kCz. Defining 

d„{t) = (d‘(r), dlit), dlit), dlit)^ (27) 

and using the expansion (25), the Dirac equation reads 

idnit) = (a • (nke^ + /?o) + l3mo)dnit) 


- w(t)qa • a(t)dn-i(t) - w{t)qa • a{ty dn+\{t) (28) 

in momentum space. Here, the window function 

fsin^^l^ if-AF^t^O, 

= if0<t<T, (29) 

I sin^ if r < r < r + Ar, 


Choosing _^^,|(r) as the initial state ait = tin = -At, the de¬ 
sired transition matrix element for transitions from the negative- 
energy state (m, T) to the positive-energy state (n, T) reads in 
terms of the Fourier space coefficients 

G( = ^lke^+pQ^'[^m(T) (33) 

and analogously for other transitions. The average number of 
created electrons with momentum nke^ + Po and spin up/down 
is given by 

A^n,yi = ^ IGCl-yy/i-myf + \GC\-)n,Vi‘m,if' • (34) 

m 

Thus, the spin-resolved average total number of created elec¬ 
trons becomes 

= y Vi.T/i ’ (35) 

n 

and the average total number of created electrons equals 

M = N^+Ni. (36) 

The Dirac equation in momentum space (28) is an infinite 
system of ordinary differential equations. In a numerical so¬ 
lution, however, we have to limit the system to a finite num¬ 
ber of modes. This is justified because the interaction term 
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FIG. 3: (color online) Time evolution of the initial momentum eigen¬ 
state Wp=o,T its spreading into neighboring positive and negative 
momentum eigenstates for cr_ = cr+ = tt/A, co - 0 All mo and ^ = 1/2. 
The vertical axis specifies the interaction time with the external field, 
and the horizontal axis specifies the momentum mode p - nke^ + 0. 
In field: The propagated state is decomposed into free momentum 
eigenstates without turning off the laser field. This is not unambiguous 
for physical interpretation. In vacuum: The external field has been 
turned off smoothly before deeomposing the propagated wave funetion, 
allowing for a physical interpretation as momentum eigenstates. 



laser halfcycles 


FIG. 4: (color online) More detailed plot of the difference between 
taking an expectation value, while the electromagnetic field is still 
acting (in field) or has been smoothly turned off (in vacuum). Both 
lines correspond to the plotted weights c„_ (38b) of Fig. 3 with mode 
number 1. Taking the expeetation values in the field leads to no unique 
physical identification and also shows faster oscillations on the order 
of the laser frequency cx). In contrast, a measurement after the field 
has been turned off yields a clear and smooth Rabi oscillation with 
respect to the interaction time with the external laser. 


w(t)qa • a(t) does not depend on the mode number n, but the 
free Dirac Hamiltonian grows asymptotically linearly with n. 
Thus, for large n transitions between different free-particle 
states are suppressed, and consequently large-^ states do not 
contribute to pair creation. The system of ordinary differential 
equations (28) is solved via a Fourier split operator method 
similar to the method presented in Ref. [44] . 


5. Numerical results 

5.1. Single-mode dynamics 


To help the reader understanding the following simulations, we 
would like to start with an example of a single time evolution. 
Both lasers are set up with linear polarization, i. e., = a+ = 

711 A, and ^ = 1/2, where ^ is determined with respect to a single 
laser field. Taking into account the electromagnetic fields of 
both lasers would yield larger values of in particular ^ = 1 in 
case of linear polarization. To calculate the expectation value 
of an electron being created with zero momentum and spin 
up, the initial state must be set to Up={)y, see Eq. (31a). Now, 
according to Eq. (21), the propagation back in time needs to be 
calculated. The duration of the external electromagnetic field 
can be varied and will be given by the number of half-cycles 
of the lasers period 


T = (number of half-cycles) x — . (37) 

OJ 

During this evolution (including the turn on [45]), the initially 
localized wave function spreads in momentum space to neigh¬ 
boring “modes” with momentum p = nke^ + /?o- Furthermore, 
the wave function will also start to occupy negative energy 
modes Let (p{t) be the propagated wave function; then 

c„+(0 = + \{up,i\ip{t ))\^, (38a) 

Cn-(t) = + \{Vp,i\(p{t))\^ (38b) 


define the weights (probability coefficients) of this wave func¬ 
tion to be in a positive (c„+) or negative (c„_) energy mode with 
momentum p = nke^ + Po. Here, the spin degree of freedom 
has already been summed out. This decomposition into free 
momentum eigenstates is not suited for physical interpretation 
[21] if the external field is still acting, because the free mo¬ 
mentum eigenstates do not represent eigenstates of the Dirac 
Hamiltonian in the presence of an external electromagnetic 
field. Hence, for a unique physical interpretation, the external 
field is turned off smoothly, and then the decomposition into 
free momentum eigenstates is carried out. Both cases, with 
and without turn off, are shown in Fig. 3. They look quite the 
same for this set of parameters, but indeed they are different. 
This difference is visualized in more in detail in Fig. 4 for the 
weight Cn- with mode number n = \. The measurement in the 
field contains oscillation noise with frequencies on the order of 
the laser frequency. These unphysical oscillations vanish, if the 
measurement is performed in a properly switched-off field. For 
the measurement in vacuum, all plots in Figs. 3 and 4 show a 
slow oscillation with respect to the duration of the electromag¬ 
netic field, corresponding to a Rabi oscillation of A/o,t, which 
is the expectation value for the creation of an electron with zero 
momentum and spin up, given by the sum over all weights c„_ 
in the negative-energy subspace. Similar Rabi oscillations have 
been predicted for the relativistic Kapitza-Dirac effect [46, 47]. 








6 


1.0 

0.5 

0.0 


(jj - OAllmoc^lH a± = 1.0 X n/i 




0 50 100 150 200 250 300 350 400 

laser halfcycles 


FIG. 5: (color online) Rabi oscillations for the average occupation 
number A/o,m electrons with final momentum p - 0. Three dif¬ 
ferent parameter sets for the polarization parameters cr_ and a+ have 
been chosen, while ^ - 1/2 is kept for all. Top: The clearly visible 
resonant Rabi oscillation is independent of the final spin orientation 
due to the linear polarization. Middle: Only the electron with spin up 
is at resonance for this specific elliptical polarization set, while the 
electron with spin down is completely off resonant. Bottom: Both spin 
orientations are slightly off resonant with almost the same oseillation 
amplitude but different Rabi frequeneies. 


5.2. Rabi oscillations 

It was mentioned in the previous subsection that the expectation 
values for creating electrons and positrons with specific mo¬ 
mentum and spin (given by the occupation number (34)) will 
exhibit Rabi-type oscillations with respect to the interaction 
time in the external field. Choosing a particular set of values 
for the polarization parameters a+ and cr_, the field strength 
(determined via ^), and the final momenta and spin (/;, 'f/^) for 
the created particles of interest, pair creation becomes resonant 
for specific frequencies of the laser field. Three distinctive 
parameter sets are shown in Fig. 5. For all plots, the parameter 
^ equals 1/2, and the final momenta of interest for the created 
electron is set to = 0. The upper plot shows that the resonance 
of the average number of created electrons is independent on 
their final spin in a laser setup with linear polarization. For 
a+ = a- = 0.6 X 7r/4, however, a spin effect can be observed 
as illustrated in the middle part of Fig. 5. The created electron 
with spin up is at resonance, while the electron with spin down 
is clearly off resonant and does not show any noticeable Rabi 
oscillation. The bottom plot shows that for an elliptical polar¬ 
ization (a+ = a- = 0.8 X 7r/4, not far from linear polarization) 
both Rabi oscillations of the two possible spin orientations are 
slightly off resonant and almost equal in their amplitude, but 
their Rabi frequencies differ. 


5.3. Resonance spectra 

In this and the following two sections, we will consider parti¬ 
cles having zero momentum only, while Sec. 5.6 will allow for 
the case of nonzero momentum. Having examined the Rabi 



FIG. 6: (color online) Pair production resonances for linear polar¬ 
ization a- = a+ = njA and ^ = 1/2. For each laser frequency du, a 
simulation of the type in Sec. 5.2 is performed, and the maximum 
value of A/o,T/i is taken as the ordinate value for this frequency. The 
spin of the produced electron needs not be distinguished for electrons 
with zero momentum due to the linear polarization. Four distinct 
resonance peaks can be seen with increasing resonance width as the 
laser frequency increases. The smaller peaks have a very low Rabi 
frequency and do not reach a full Rabi cycle as the maximum interac¬ 
tion time of this simulation is limited to 401 half-cycles of the laser 
field. For longer interaction times, these resonance peaks will grow, 
and new ones may emerge. 
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FIG. 7: (color online) Resonance at du ^ 0.472 mo for linearly po¬ 
larized waves with a+ = a- = njA and ^ = 1/2 in detail. Left: The 
pair production amplitude resonance including a Lorentz fit (39) for 
the resonance peak. Right: The Rabi oseillation frequeney Q and its 
dependence on the detuning, including a fit of the form (40). Both fits 
agree perfectly well. The numerical fit values for the Lorentz fit are 
du^ = 0.47174(1) mo and w = 0.001780(3) mo, and a fit to (40) yields 
= 0.002012(4) mo, rj = 2.26639(4), and du^ = 0.47175(1) mo. 


oscillations for specific laser frequencies, we will analyze the 
pair-creation spectra now. Each data point of a pair-creation 
spectrum represents the maximum value of A/o,t/I for a spe¬ 
cific set of laser parameters. Figure 6 shows the maximum 
value of JVo,yi as a function of the laser frequency du while 
keeping all other parameters fixed. Here the polarization is 
linear, i. e., dr_ = a+ = njA, and ^ = 1/2. Note that in case of 
linear polarization pair-creation probabilities do not depend 
on the spin orientation, and therefore A/o,t = A/o,|. Four com¬ 
plete resonances can be inferred from Fig. 6. Their resonance 
width increases with increasing frequency. Furthermore, three 
smaller resonances appear that do not reach the full resonance 
maximum of unity. The reason for this is that their Rabi oscilla¬ 
tion frequencies are very low and the top of the Rabi cycle has 
not been reached before the simulation ends as the maximum 
interaction time with the external field has been limited to 401 
half-cycles of the external laser field for this simulation. 

A close-up sweep of the resonance at du ~ 0.472 mo has been 
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FIG. 8: (color online) Polarization and laser frequency dependence of the electron creation probability with zero final momentum and specific 
spin orientation with respect to the positive z direction. Both lasers field strengths’ are determined by their frequency cx) and ^ = 1/2. The 
color encodes the maximum electron creation probability measured for a time of maximum 400 half-cycles of the laser fields period. Various 
resonances can be seen, which change their resonance frequency ojr with varying ellipticity a+. Left two plots: Both plots correspond to the 
case of corotating laser fields. A slice along the line of linear polarization a+ - n/A would correspond to the spectrum Fig. 6. Except for 
Qf+ = njA, which corresponds to linear polarization, electron creation probabilities depend on the spin orientation of the created electron. Note 
that pair creation is completely suppressed for circular polarized light (a+ = 0 or cr+ = njl). Right plot: This represents the resonance spectra for 
antirotating fields, where the electron creation probabilities do not depend on the spin orientation. 


done to show the typical Rabi characteristics. One characteris¬ 
tic, the Lorentz type resonance shape, is plotted on the left of 
Fig. 7 with a fit to the Lorentz function 

maxA/o,t/i =-^— 

l + (^ 

The fit parameter ojr corresponds to the resonance frequency, 
and w gives to the full width at half maximum. Furthermore, 
the Rabi oscillation frequency O can be extracted from the 
numerical data for each laser frequency oj. The right part of 
Fig. 7 shows the Rabi frequency O together with a fit to the 
function 


n= (40) 

where CIr denotes the Rabi frequency at resonance and 77 is a 
fit coefficient depending on the resonance. It can be seen that 
both fits agree perfectly well with the numerical data and hence 
justify the categorization of pair production in this parameter 
regime as a Rabi process. 

Extrapolating these results into the regime of frequencies 
lower than shown in Fig. 6 , the resonance peaks become very 
dense with respect to the frequency axis. This is due to the fact 
that the positions of neighboring resonances, which differ in 
the number of absorbed photons by 1 , will differ only by the 
small photon energy corresponding to the low laser frequen¬ 
cies. The Rabi frequencies of the resonances will decrease for 
the increasing number of photons needed to create a pair; see 
Ref. [21]. Hence, only the very beginning of the Rabi oscil¬ 
lations and thus a quadratic time dependence of the average 
occupation number for a created pair with sharp momentum 
at a specific resonance will be observable. Incorporating the 
finite resolution of measuring the momenta of the produced 
pairs and the corresponding integration in momentum space 
(analogous to Fermi’s golden rule), we expect that this leads to a 


linear time dependence, i. e., a rate for pair production, which is 
characteristic for low laser frequencies in the tunneling regime. 


5.4. Polarization dependence of the 
resonance spectra 

Going one step further, the resonance structure in Fig. 6 can also 
be calculated for different polarizations than a+ = a- = 7 r/ 4 . 
Performing a two-dimensional sweep over the laser frequency 
OJ and the polarization parameters a+ and a- leads to the two- 
dimensional color plots shown in Fig. 8 . The polarization 
parameters a+ and a- are not varied independently; instead 
they either fulfill a+ = a- (corotating fields) or a- = n/l - a+ 
(antirotating fields). 

Figure 8 shows the spin resolved electron creation proba¬ 
bility for electrons with zero final momentum as a function of 
the laser frequency oj and the polarization. The two plots on 
the left correspond to the spin-up and spin-down spectra for 
corotating fields, while the single plot on the right presents the 
spectrum for antirotating fields. Our numerical results show 
that the spectra depend on the spin orientation for corotating and 
are identical for both spin orientations in case of anti-rotating 
fields. Because spin-down and spin-up spectra are different 
for corotating fields, see left part of Fig. 8 , it is possible to 
have for some fixed oj and a+ zero electron creation probability 
for a specific spin orientation and nonzero probability for the 
opposite orientation. For antirotating fields, however, the elec¬ 
tron creation probability is completely symmetric with respect 
to the spin orientation of the created electron as shown in the 
right part of Fig. 8 . Note that pair creation is suppressed for 
a+ = 0 and a+ = njl, i. e., circular polarization, in the case of 
corotating fields. For antirotating fields, however, resonances 
exist also for a+ = 0 and a+ = njl. 

A possible interpretation of the different behavior of corotat- 
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ing and antirotating fields is that the standing electromagnetic 
wave carries nonzero spin density for corotating fields but not 
for antirotating fields. The orientation of the photonic spin den¬ 
sity breaks the symmetry with respect to reflection along the 
propagation axis of the laser fields; see Fig. 2. Consequently, 
pair creation probabilities depend on the spin orientation of the 
created electrons. Note that for a+ = njA (linear polarization) 
the photonic spin density vanishes also for corotating fields. 
Consequently, there is no dependence on the electron spin; 
see left part of Fig. 8. Pair creation probability becomes spin 
dependent, however, if a+ njA. 

The parameter a+ is restricted to the range from 0 to njl 
in Fig. 8. One can show by the symmetry properties of the 
standing electromagnetic wave’s vector potential (or its elec¬ 
tromagnetic held) that for a+ = a- the spectra in Fig. 8 must 
be periodic in a+ with period tt. The part of the spectrum for 
7112 < a+ <71, which is not shown in Fig. 8, can be obtained 
by mirroring the spectrum for 0 < a+ < 7r/2 at a+ = Till. Fur¬ 
thermore, the spectrum is periodic in a+ with period 7r/2 for 
antirotating fields. 

Figure 8 reveals an explicit spin dependence of the reso¬ 
nance structures in case of corotating fields. These differences 
become more pronounced the more the fields’ ellipticity differs 
form linear polarization at a+ = TijA. The spectra for spin up, 
if mirrored along the line a+ = TijA, are equal to the spectra 
for spin down and vice versa. This can be understood by sym¬ 
metry, because flipping the spin of the final particle leads to 
the same result as flipping the sense of rotation for the lasers’ 
polarization, i. e., from right circulating polarization to left 
circulating polarization. Looking at the resonance between 
0.8 mo < u; < 0.9 mo, for example, the resonance frequency ojr 
decreases toward circular right polarization a+ Tijl for spin 
up, while it increases for spin down. 

Note that this shift of the resonance frequency is always 
opposite for electrons with spin up and with spin down. This 
might be understood by the coupling between the particle’s 
spin and the lasers’ magnetic field, as seen in the rest frame 
of the particle. With the electric and magnetic field compo¬ 
nents E and B in the laboratory frame, the magnetic field in 
the rest frame is in leading order B' ~ B - v x E. Because 
B is always perpendicular to the z axis, the Zeeman energy 
-ql{2mo) cr ' B' = ql(2mo) cr • {vxE) results, where v denotes 
the particle’s velocity in the laboratory frame. For corotating 
laser fields with elliptical polarization, the electric field rotates 
in the x-y plane and therefore leads to a two-dimensional tra¬ 
jectory of the particle in the x-y plane. Because E also lies in 
the x-y plane, the magnetic field vxE points in the z direction 
yielding a modification of the particle’s energy and in this way 
modifying the resonance frequency ojr. The sign of this energy 
shift depends on the particle’s spin in the z direction, which 
explains the opposite shift of the resonance frequency for dif¬ 
ferent spin orientations. For antirotating fields, in contrast, the 
electric field does not rotate in the x-y plane. It oscillates with 
a fixed direction in the x-y plane. Consequently, the particle’s 
trajectory projected onto the x-y plane will always be parallel 
to the electric field direction. Thus, the Zeeman energy is zero 
for an electron with spin in z direction, and consequently there 
is no spin-dependent shift of the resonance frequency. 


The shift of the Zeeman energy, which we find in the case of 
corotating fields, is also given by the nonrelativistic corrections 
to the Pauli equation. These correction can be obtained by an 
expansion of the Dirac equation via a Foldy-Wouthuysen trans¬ 
formation and contains the term l/(4mQ) (r ' (E x (-iV - qA)). 
For zero canonical momentum, this is up to a factor of 2 equiv¬ 
alent to the shift of the Zeeman energy discussed above. Note 
that the factor E x A can be attributed to the photonic spin 
density of the electromagnetic field [34], which is zero for 
antirotating fields. 

5.5. Comparison to the dipole 
approximation 

When pair creation in standing electromagnetic waves is stud¬ 
ied, neglecting the magnetic field component and the spatial 
variation of the electric field is a commonly applied approxi¬ 
mation [31, 48, 49]. Therefore, we will analyze the effect of 
this approximation in this subsection. Neglecting the spatial 
variation of the electromagnetic field given by the potential (4) 
yields the dipole approximation vector potential 

A+{t) = Re (cos a+ \l) + sin a+ e^*^- |r)^ e"^^^j. (41) 

This dipole approximation of the electromagnetic fields is em¬ 
ployed in Fig. 9 for comparison with Fig. 8, where the full elec¬ 
tromagnetic field was taken into account. Also in the framework 
of the dipole approximation, Rabi oscillations of the electron 
creation probability and resonances can be found. The dipole 
approximation gives, however, results, which differ quantita¬ 
tively as well as qualitatively from the case of pair creation in 
electromagnetic fields. For example, the number of resonances 
is reduced in case of the dipole approximation. In Fig. 9 the 
parameter a+ ranges from 0 to tt in order to point at differences 
to the simulation with the full monochromatic field and to show 
some peculiarities of the turn on and off of the electromagnetic 
field. 

In case of corotating fields, the two plots on the left of Fig. 9 
are not perfectly symmetric with respect to the line a+ = 7r/2, 
especially for 0.9 mo < oj < 1.2 rriQ. This is due to the turn on 
and turn off as the following argument illustrates. Comparing, 
for example, a+ = a- = TijA with a+ = a- = 37il4, both laser 
fields are linearly polarized with the electric field component 
pointing in the v direction and in the y direction, respectively. 
Due to rotational invariance around the z axis, both parameters 
should yield the same results. But owing to the parametrization 
and analytical structure of the vector potential (41), the electric 
field amplitude is at maximum for a+ = a- = TijA at r = 0, 
while it is zero for a+ = a- = 37rl4 slI t = 0 maximal. This 
phase shift leads in combination with the window function (29) 
to an asymmetry, which has no effect in regimes with Rabi 
frequencies much slower than the laser frequency, because the 
pair production effect builds up over many laser cycles and is 
less sensitive to the initial phase of the laser field. For Rabi 
frequencies on the order of the laser frequency as observed for 
oj> mo this difference, however, becomes noticeable. 
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FIG. 9: (color online) Polarization and laser frequency dependence of the electron creation probability with zero final momentum and specific 
spin orientation with respect to the positive z direction when the dipole approximation is applied. All parameters as in Fig. 8. 
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FIG. 10 : (color online) Polarization dependence of the pair creation probability, where the electron has different final momenta po = (0, 0,Pz). 
The lasers’ frequencies are set to u; = 0.8 mo, and their strength is implicitly determined by ^ = 1/2. From all four plots can be inferred that, for 
this fixed frequency and strength, final pairs with distinguished momenta will be produced. Left two plots: These two plots correspond to the 
corotating field setup and differ by the spin orientation of the final electron. The vertical cut along Pz = 0 corresponds to a vertical cut along 
oj = 0.8 mo in Fig. 8. Both plots are symmetric around p^ = 0. This is because inverting p^ would also invert the final helicity of the created 
particles, but due to the fact that for corotating fields both fields’ helicities are opposite to each other, the systems stays invariant. Right two plots: 
These two plots correspond to the antirotating field setup. In contrast to the case of corotating fields, the spectra are not invariant under inversion 
of the Pz momentum and consequently the helicity of the final particle, because both laser fields have the same helicity for this case. But if, in 
addition to the momentum inversion, also the spin is inverted, the helicity of the final particle does not change, and the spectrum stays the same. 


For the case of antirotating electromagnetic fields, we found 
in Sec. 5.4 that the probability to create an electron with zero 
momentum is periodic in the field parameter a+ with period 
7r/2. Applying the dipole approximation, however, the spectrum 
becomes periodic in a+ with period tt as shown in the right 
part of Fig. 9. This is a direct consequence of the symmetry 
properties of the vector potentials in (41). For a+ = njA there 
are broad resonances, whereas pair creation is inhibited for 
a+ = ?>7 iIA. Both parameters correspond to electric fields with 
linear polarization. For a+ = njA the vector potentials A+{t) 
and A-{t) given in (41) superimpose constructively, while both 
potentials cancel for a+ = 37r/4. Note that this cancellation is a 
peculiarity of the dipole approximation. Two counterpropagat- 
ing electromagnetic waves with temporal and spatial variation 
as specified in (4) cannot interfere such that both annihilate 
each other completely. 

Notably, the simulation for the dipole approximation also 
yields a strong dependence of the pair creation on the spin of 


the created electron and the polarization of the electric field for 
corotating fields while not for antirotating fields. Comparing 
the magnitude of the shift of the resonance frequency ojr with 
the one in Fig. 8 shows that this effect is even stronger in the 
dipole approximation with a rotating electric field only. As 
mentioned in the previous subsection, this shift might be due 
to the coupling of the particles spin and the magnetic field, as 
seen in the rest frame of the particle. Alternatively, the spin 
dependence could also be inferred from symmetry. In case of 
corotating electric fields, the sense of rotation changes under 
reflection of the z axis, and thus the pair creation probability 
might become spin dependent. In contrast, two antirotating 
electric fields form an electric field with a constant direction, 
and therefore the field configuration is symmetric under reflec¬ 
tion of the z axis. Consequently, also pair creation is symmetric 
under reflection of the final particle’s spin, assuming zero final 
momentum. 

Comparing the various resonances in Figs. 8 and 9 in terms 
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of the polarization dependent shift of the resonance frequencies 
shows that the frequency shift becomes smaller at lower laser 
frequencies. The reason for this might be a combination of the 
decreased amount of photon energy and an averaging effect 
due to the increased number of photons. The polarization 
dependence of the resonance frequencies tends to be larger for 
the dipole approximation than for the case of electromagnetic 
fields, which may be attributed to the fact that the electric field is 
maximal at each point in space within the dipole approximation 
and to the absence of the magnetic field in the laboratory frame. 


5.6. Final momentum dependence of the 
resonance spectra 

In the previous subsections, we investigated the creation of 
particles with zero final momentum. Looking, for example, at 
Fig. 8, it shows that at Cc) ~ 0.8 mo no particles with zero final 
momentum can be created for the given setup. Particles with 
nonvanishing momentum, however, might be created. A sweep 
over the final momentum allows us to study the momentum 
dependence of the pair production probability. Due to the fact, 
that a three-dimensional sweep over all final momenta is nu¬ 
merically quite expensive and also hard to visualize, we vary 
only the momentum directed along the z axis, i. e., the lasers’ 
propagation axis. The perpendicular momenta are set to zero. 
Furthermore, we chose a fixed laser frequency of u; = 0.8 mo 
and additionally scanned over the ellipticity for corotating and 
antirotating field configurations as done in Figs. 8 and 9. The 
numerical result is shown in Fig. 10. The two plots on the left 
correspond to corotating fields, while the two on the right cor¬ 
respond to antirotating fields. In agreement with our previous 
calculations. Fig. 10 indicates that no particles are created with 
zero final momentum for the fixed frequency of cl> = 0.8 mo- 
Pair creation is possible for some nonvanishing momenta as 
the resonances in Fig. 10 indicate. A notable peculiarity for the 
case of corotating fields can be seen in the left part of Fig. 10. 
Tracing the two broader resonances close to = 0 for the spin 
up case from top to bottom, they approach each other but do 
not join each other at = 0 for a+ ~ 7r/8. Thus, pair creation 
is suppressed at = 0 exactly consistent with Fig. 8. 

These resonances in Fig. 10 exhibit distinctive symmetries. 
For corotating fields, both plots for spin up and spin down are 
symmetric around the line Pz = 0. This symmetry is related 
to the fact that inverting the z momentum preserves the rela¬ 
tive orientation of the helicities of the created electron and the 
two counterpropagating waves. Inverting the final momentum 
inverts the helicity of the electron. Because both laser fields’ 
helicities are opposite to each other for corotating fields, in¬ 
verting the z momentum or equally the helicity of the created 
electron does not change the system. The electron’s helicity 
is parallel to one wave and antiparallel to the other as before 
inversion of the z momentum. 

Both helicities of the laser fields point into one preferred di¬ 
rection for antirotating fields. Thus, inverting the z momentum 
does not leave the system invariant, and electrons are created 
preferably with helicity parallel or antiparallel to the helicity of 


the electromagnetic fields. To preserve the relative orientation 
of the helicities of the created electron and the two counterprop¬ 
agating waves, one has to invert the z momentum and the waves’ 
helicity, which corresponds to the combined transformation 
Pz Pz ^ thus the spectra in the right 

part of Fig. 10 are invariant under this transformation. 

Both the corotating as well as the antirotating case have in 
common that mirroring the spectrum for one spin orientation 
around the a+ = njA line yields the spectrum for the other spin 
orientation. This symmetry holds for the spectra in Fig. 8 as 
well as in Fig. 10. Again, this is a consequence of the system’s 
invariance under flipping the created electron’s spin/helicity 
and the field’s spin/helicity. 

Let us compare the spectra for corotating fields with the ones 
for antirotating fields in Fig. 10 at fixed values of the ellipticity 
parameter a+. In the case of corotating fields, electrons with 
spin up or down are created with higher probability depend¬ 
ing on the ellipticity parameter a+ \ see the left part of Fig. 10. 
Because of the symmetry properties described in detail above, 
the electrons are created with equal probability having positive 
or negative helicity; this means created electrons propagate in 
the positive or negative z direction with equal probability. This 
changes if we consider antirotating fields; see the right part of 
Fig. 10. Here electron creation is spin symmetric. This means 
spin-up and spin-down electrons are created with equal prob¬ 
ability. However, electrons with positive or negative helicity 
are created with higher probability depending on the elliptic¬ 
ity parameter a+. For an ellipticity of a+ ~ 7r/8, for example, 
electrons are created preferably with spin up, which propagate 
only along the positive z direction, and with spin down, which 
propagate only along the negative z direction. 


6. Conclusions 

We studied pair creation in strong electromagnetic fields in 
the nonperturbative regime, i. e., ^ ~ 1. For this purpose, we 
derived a numerical framework for solving the time-dependent 
Dirac equation for an electron in two counterpropagating 
monochromatic laser fields with equal frequency. This setup is 
quite generic, because also other setups with two lasers having 
unequal frequencies or a colliding angle different than 180° 
could be incorporated by a proper Lorentz transformation. Fur¬ 
thermore, the polarization and the field strength of each in¬ 
dividual laser field could be chosen arbitrarily. With this we 
overcame the so often applied dipole approximation for the laser 
fields and included the electric field’s space dependence as well 
as the magnetic field’s influence on pair production, leading 
to true three-dimensional dynamics. For numerical feasibility 
and in order to keep the space of possible parameters small, 
high laser frequencies and equal intensity and polarization of 
the counterpropagating fields were assumed. The numerical 
method, however, is also capable of treating setups with un¬ 
equal polarizations and intensities. It can be easily generalized 
to multicolor setups, provided that all frequencies are low-order 
harmonics of some base frequency. 

It was shown that switching off the external laser fields is 
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essential for a proper determination of the pair creation proba¬ 
bilities. These probabilities exhibit the characteristic features of 
Rabi oscillations with the corresponding resonance frequencies 
and Rabi frequencies. The resonance structures and their de¬ 
pendencies on the final quantum numbers of the produced pairs 
as well as on the polarization of the laser field were investigated 
numerically in more detail. 

A quite remarkable result is the observation that, by tuning 
the polarization of the two laser beams, spin polarized pairs can 
be created. By using corotating laser fields, for example, the 


spins of the electrons in the produced pairs may be aligned all 
the same along the propagation axis of the lasers. Antirotating 
fields, in contrast, may yield electrons with spin-up propagating 
forward and electrons with spin down propagating backward 
with respect to the lasers’ propagation axis. 
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